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ABSTRACT 



Unpolarized Gowdy models are inhomogeneous cosmological models that depend on 
time and one spatial variable with complicated nonlinear equations of motion. There 
are two topologies associated with these models, T 3 (three-torus) and 5" 1 x 5" 2 . The 
T 3 models have been used for numerical studies because it seems to be difficult to find 
analytic solutions to their nonlinear Einstein equations. The 5" 1 x 5" 2 models have even 
more complicated equations, but at least one family of analytic solutions can be given as 
a reinterpretation of known solutions. Various properties of this family of solutions are 
studied. 



I. Introduction 



In 1971 Gowdy proposed a set of inhomogeneous cosmological metrics with the idea 
of constructing vacuum models in which inhomogeneity appeared in a simple way [1]. 
He wanted his models to have compact topology and two commuting Killing vectors con- 
structed so that the model depended only on time and one of the spatial variables. Models 
of this sort have the metric structure of axisymmetric or cylindrically symmetric static or 
stationary metrics with the radial coordinate replaced by t and the metric independent 
of the new "radial" coordinate instead of being independent of time. Such metrics share 
many of the properties of axisymmetric metrics and a number of well-known results for 
these metrics can be applied to the Gowdy models almost without change. 

There are two topologies for the t = constant three surfaces of these metrics that were 
proposed by Gowdy. One was a T 3 (three-torus) topology and another was an 5" 1 x S 2 
topology. For some reason the T 3 topology metrics have been used extensively for a number 
of studies of different types, while the models with S 1 x S 2 topology have been almost 
forgotten. 

In the notation of Berger and Garfinkle [2] the original T 3 Gowdy model takes the 

form 

ds 2 = e- A /V/ 2 (-e- 2T dr 2 + d9 2 ) + e~ T [e p da 2 + e~ p dd% (1.1) 



where A and P are functions of 9 and r. The model is compactified by requiring < 
6, a, 5 < 2/T, and the metric functions must be periodic in 9. It is often stated that 
this form of the metric can be achieved by coordinate transformations alone, but (see for 
example, [3]) the form of the coefficient of dd 2 comes from the fact that the determinant of 
the two-metric g 2 formed by the da 2 and d5 2 terms, which in this case is e~ 2r , is a solution 
of a certain combination of the Einstein equations. If we assume that g 2 is g 2 (T,9), and 
define T = e _T , then 
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This expression is valid for any two- metric analogous to the da 2 , dS 2 part of (1.1) of the 
form A(9, r)da 2 + 2B(9, r)dad5 + 0(0, r)d5 2 , not just the special case given in (1.1), a fact 



that we will use below. In the case of (1.1) <72 = T 



-It 



is a solution of (1.2), which 



allows us to write the dd 2 term of (1.1) in the form given above. 

The Einstein equations for (1.1) become a linear equation for P and a set of first-order 
equations for A which can be integrated to give A given any solution to the equation for 
P. These Einstein equations are easily solved in terms of ordinary Bessel functions, and 
several particular solutions were given by Gowdy [1]. 

Gowdy [4] also proposed "unpolarized" models (the name comes from the fact that for 
small values of P the metric [1.1] becomes that of polarized linearized gravitational waves, 
and the unpolarized models become that of unpolarized waves), which, in the notation of 
[2] become, for the T 3 models, 



ds 2 = e- x / 2 e T / 2 (-e- 2T dr 2 + d9 2 ) 



+e~ T (e p da 2 + 2e p Qdad5 + [e~ y + e p Q 2 ]d5 2 ), 



(1.3) 



where the form of the d5 2 term comes from taking g 2 = e 2r as above, and A, P, Q 
are functions of 9 and r. The Einstein equations for (1.3) are a set of coupled nonlinear 



1 



equations for P and Q, and, again, first-order equations for A that can be integrated once 
P and Q are known. 

The models (1.3) have been used for numerical studies [5], since it is usually assumed 
that there is little chance of finding analytic solutions for such a complicated system of 
equations. One aim of these studies has been to show that, in general, the solutions to 
these models become asymptotically velocity term dominated (AVTD) as r — > +00, that 
is, near the singularity of the model. The idea of AVTD models, originally due to Belinskii, 
Khalatnikov and Lifshitz [6] (BKL), is that near a singularity in the ADM action, 

I = 7^— / [?r% - NH ± - NiH^dtd^x, (1.4) 



re = -2^., (i.56) 

where 3 R is the curvature on t = const, surfaces and | is a covariant derivative on these 
surfaces, the W = terms will imply that the curvature ("potential") term in the H± 
can, at best, be ignored completely, or at worst become algebraic in g^, so that the 
ADM equations can be solved as if the metric were homogeneous (with the completely 
independent functions of time at each point x). The conjecture of BKL is that this behavior 
is generic, an idea which is still being debated. Numerical studies of the models (1.3) (see 
[2] and references therein) have shown that these models do indeed become AVTD except 
at a set of isolated points (9 = where there are "spikes" in the solution. The origin of 
these spikes is studied in Ref. [2]. 

While the T 3 Gowdy models have been used extensively in classical and minisuper- 
space quantum gravity [7] , the S 1 x S 2 topology seems to have been almost ignored. The 
polarized form of these models can be written, in the notation of [2], as 

ds 2 = e- x ' 2 e T ' 2 {-e- 2T dr 2 + d9 2 ) 

+ sin(e- T )[e p cM 2 + e" p sin 2 6># 2 ], (1.6) 

where gi = sin 2 T sin 2 9 is a solution of (1.2). These models have Einstein equations that 
consist, as above, of a linear equation for P and first-order equations for A that can be 
directly integrated. The linear equation can be solved by separation of variables, and 
several particular solutions were given by Gowdy [1]. 
The unpolarized form of these models, 

ds 2 = e- x ' 2 e T / 2 {-e- 2T dr 2 + d9 2 ) 

+ sin(e- r ) [e p d5 2 + 2e p Qd8d(j) + (e p Q 2 + e~ p sin 2 6>)# 2 ], (1.7) 

(where, once again, gi = sin 2 e~ T sin 2 9 gives the form of the d(p 2 term), has Einstein 
equations similar to those of (1.3), though quite a bit more complicated. These models 
seem to have been ignored, and there are, as far as we know, no numerical or analytic 
solutions given in the literature. 

One could argue that the fact that the Einstein equations for (1.7) are even more 
complicated than those for (1.3) would mean that only numerical solutions could be found. 
However, we will show that at least one family of particular analytic solutions is readily 
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available as a reinterpretation of the portion of the Kerr metric [7] between the inner, 

r_ = M — V M 2 — a 2 , and the outer, r+ = M + \J M 2 — a 2 horizons (where the "radial" 
coordinate is timelike). We will study some of the properties of this family of solutions 
and compare some of the results with those of numerical solutions of (1.3). One important 
point is that these solutions are not AVTD. Of course, since such particular solutions 
form a set of measure zero in the space of all solutions, it is not clear what the generic 
behavior of these models near the singularities is. Another point is that the "singularities" 
of these models are not true curvature singularities, since they represent points where the 
solution crosses a lightlike surface into regions where the metric is the usual stationary 
axisymmetric black hole metric. We will discuss possible generalizations briefly. 

The paper is organized as follows: in Sec. II we will give the Einstein equations for all 
the metrics listed above, and show that the Kerr metric between its horizons is a family of 
exact solutions for the unpolarized S 1 x S 2 model. Sec. Ill discusses the properties of the 
solutions (including possible AVTD behavior). Sec. IV consists of conclusions and some 
suggestions for extending the results of this article. 
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II. Exact Solutions 

Equations of Motion 

The equations of motion for the polarized Gowdy models were given in the original 
works of Gowdy [1,4]. In our parametrization the equations for the T 3 cosmologies become, 
as we have mentioned above, a linear equation for P, 

P, TT -e- 2T P,gg = 0, (2.1) 

and a set of first-order equations for A, 

X, T = P 2 T + e- 2T Pl (2.2a) 

A «, = 2P,gP jT , (2.26) 

which, in principle can be integrated directly to give A if they obey an integrability con- 
dition which is just (2.1). Eq. (2.1) is exactly soluble in terms of trigonometric functions 
of 9 and ordinary Bessel functions of r, and there are even tabulated integrals that give A 
explicitly in certain cases [8] . Some of these solutions were given by Gowdy in his original 
papers. 

The equations for models with S 1 x S 2 topology are only slightly more complicated, 
again we have a linear equation for P, 

P,rr ~ ^— ksinflPfl),* - e~ 2r - [e~ T cot( e - T ) - 1]P T = 0, (2.3) 
sin u 

this equation can be separated in terms of Legendre polynomials and a function of r. Some 
exact solutions to this equation were given by Gowdy [1] (in another coordinate system). 
The equations for A become 

cot(e- T )A, e - 2e T (P T P fl ) + cot 6[-e T \ iT + 2e T P T +e T +2 cot(e" T )] = 0, (2.4a) 

cot(e-)(A, T -l)-e-[(P r ) 2 + e- 2 -(P 9 ) 2 ] 
+e~ T [cot 2 (e" r ) + 4] + e" T (- cot 6\, + 2 cot 0P e ) = 0. (2.46) 

The unpolarized models have somewhat more complicated equations. The equations 
for P become nonlinear because nonlinear terms in Q are added, and an equation for 
Q appears which has nonlinear terms in both Q and P. The T 3 model has been used 
extensively for numerical studies of such concepts as AVTD models [2, 5]. Its equations 
for P and Q are 

Ptt _ e - 2T P,ee - e 2P (Q 2 T - e~ 2T Q 2 e ) = 0, (2.5a) 

Q,tt - e- 2T Q, 9e + 2(P T Q, T - e- 2T P, e Q, e ) = 0, (2.56) 
and once these are solved one can intregrate the equations for A, 

A, T - [P 2 + e- 2 -P 2 + e 2P (Q 2 T + e^Q 2 ,)] = 0, (2.6a) 

X,e - 2(P e P T + e 2P Q, Q, T ) = 0. (2.66) 
4 



One of the reasons usually given for solving these equations numerically is that there is little 
chance of finding analytic solutions to such a complicated system of nonlinear equations. 

The unpolarized S 1 x S 2 models have even more complicated dynamical equations for 
P and Q, 

P,rr ~ e- 2T ( " n *y ,g - e" 2 - - -^-A(Q,r) 2 ~ (Q ,ef} - [e T cot(0 - 1]P, T = 0, 
sin a sin 

(2.7a) 

Q,rT-e- 2T Q,ee+e- 2T cot6Q, e + 2(P,TQ,T-e- 2T P, e Q, e )-[e T cot(e- T )-l]Q, T = 0. (2.76) 

There are two coupled first-order equations for A (which, of course, can be reduced to 
separate equations for A )T and A^, but give equations that are more complicated than the 
originals). They are 

cot(e- T )A e - 2e T (P T P + e 2P ^^) 

+ cot #[-e T A, r + 2e T P T + e T + 2 cot(e _T )] = 0, (2.8a) 

ip 

oot(e- T )(A, T - 1) - e T [(P T ) 2 + e^P,) 2 ] - e^*[(Q, T ) 2 + e~ 2 ^Q, e ) 2 ] 

sin 6/ 

+e- T [cot 2 (e" T ) +4] + e" T (- cot ^A, e + 2cot^P fl ) = 0. (2.86) 

One could expect that it might be difficult to find exact solutions to this system of equa- 
tions, but one is at hand without any work. 

Black Holes inside their Horizons. 

Kantowski and Sachs [9] studied cosmological models with four-dimensional groups of 
motion which had three-dimensional subgroups, and one special case was a cosmological 
model of the form 

ds 2 = -N 2 (t)dt 2 + e 2 ^^dr 2 + e -2V3/3(t) e -2V3n (t)(d(9 2 + gin 2 ^ 2) _ (2 9) 
They gave a one parameter family of solutions of the form 

t 



N 2 = - 1 , e 2 ^= --1, e -2V3(3(t) e -2V3a(t) = ^ (2 1Q) 

"7" 1 ^ 



They note that this family is nothing more than the Schwarzschild model inside the horizon. 
That is, if we write the usual Schwarzschild metric in Schwarzschild coordinates, we have 



ds 2 = -(1 - — )dt 2 + 1 2m dr 2 + r 2 (d6 2 + sin 2 Ochp 2 ), (2.11) 



and for r < 2m we can make the coordinate transformation t «-> r and we find (2.9) with 
the functions given by (2.10). 

One advantage of looking at the Kantowski-Sachs model as Schwarzschild is that we 
can see the character of the singular points of the metric. The singularity at t = is 
the true singular point of Schwarzschild (r = 0), where the curvature is infinite, but the 
singularity at t = a is just a lightlike surface where the curvature is regular and we pass 
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from the Kantowski-Sachs region to the Schwarzschild region. Something similar happens 
with the vacuum Taub model [10], where there are two singular points, but both of them 
are lightlike surfaces where we pass into the NUT [11, 12] region. 

The Gowdy models have two commuting Killing vectors, and in the T 3 and S 1 x S 2 
topologies are axisymmetric. If we were able to find a static or stationary axisymmetric 
metric with a horizon of the Schwarzschild type for which the region inside of the horizon 
had the proper character, then one could find a solution to a Gowdy model in the same 
way that one finds the Kantowski-Sachs solution from Schwarzschild. One difference is 
that in the Schwarzschild case the Birkhoff theorem says that the Kantowski-Sachs model 
is the general solution to the Einstein equations for this model. In the Gowdy case we can 
probably expect no more than a particular solution. 

A metric that generalizes Schwarzschild and has two commuting Killing vectors is the 
Kerr metric. In Boyer-Lindquist coordinates the Kerr metric is 

,2 r2 ~ 2Mr + a 2 , , 2 . , , l2 sin 2 9 2 2 . , , 

ds = 5 ^ o~7T \dt — a sin 0d(p\ +^5 5 ttt ( r + a )d(p - adt\ 

r 2 + a 2 cos 2 9 1 1 r 2 + a 2 cos 2 9 LV ' 1 

> J " 2 cos2 + ( r a + „-> cos 2 )d9 2 . (2.11) 



r 2 - 2Mr + a 2 



For this metric the outer horizon, r_|_ = M + \/M 2 — a 2 is a lightlike surface where the light 
cones tip over to the point where d/dr is a timelike vector, so that a coordinate change of 

the r <-> t is possible. Unfortunately, at the inner horizon, r_ = M 2 — \/M 2 — a 2 , another 
lightlike surface, the light cones tilt back to the point where d/dr becomes spacelike again. 
In the region r_ < r < r + we can make the transformation r «-> t and the Kerr metric 
becomes a cosmological model, 

ds 2 = ~f ~f a [dr - a 2 sin 2 9d<P] 2 + . - . J (t 2 + a 2 )# - arfr] 2 
t 2 + a 2 cos 2 6» L J t 2 + a 2 cos 2 6> J 

f2 h <r; '" 0r, ^dt 2 + (t 2 + a 2 cos 2 ^ 2 . (2.12) 



2Mt -t 2 -a 2 
This metric, with the simple transformation 

t = a y\-l3 2 cos(e- T ) + 1], (2.13) 
a = M, (3 = a/M, makes (2.12) into (1.7) with 



A = t - 2 ln(a 2 [y/l - P 2 cos(e~ T ) + l] 2 + (3 2 cos 2 0), (2.14a) 

P = ln[(l - /3 2 ) sin 2 (e- r ) + [3 2 sin 2 6»] - \n[ay/l - [5 2 sin( e - T )] 

- ln([^l-P 2 cos(e- T ) + l] 2 + [3 2 cos 2 6), (2.146) 



= 2a/jsin 2 fl[ v /r^cos(e-) + l] 
V (l-/3 2 )sin 2 (e--)+/3 2 sin 2 ^ 1 ' ' 

Even though it is obvious that Eqs. (2.14) are a solution of the vacuum Einstein equations 
(since the Kerr metric is a solution), we have, using a symbolic manipulation program, in- 
serted (2.14) into (2.7) and (2.8) and showed explicitly that these expressions are solutions. 
In the next section we will discuss the properties of this family of solutions. 
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III. Properties of the Solutions 



The general behavior of this family of solutions is given in Figs. 1 and 2, which show 

P + ln[a^/l — j3 2 sin(e~ T )] and Q as functions of 9 for various values of r. The functions 
are relatively flat at r = — hi7r, which is the "big bang", then grow a spike in the region 
of 9 = it/2 which disappears as r — > +oo, so the functions become flat again at the 
"big crunch". As r goes either to — ln7r or +oo, P becomes infinite, but the combination 

P + \n[a^l — f3 2 sin(e _r )] (i.e. adding a ^-independent term) is finite for all 9 and r except 
at the extremes of 9, 9 = and 9 = n, where, for r = — ln-zr and r = +oo, it becomes — oo 
(e p = 0). 

At first glance the "spiky" behavior of these solutions seems to have something to 
do with similar spikes that occur in numerical solutions. While we will show that there 
is a remote connection between the two situations, the spikes in our solution show up at 
different parts of the evolution from the numerical spikes. 

The simplest manner of showing the difference between the numerical and analytic 
solutions is to study the Hamiltonian formulation of the Einstein equations for the S 1 x 5" 2 
topology, which is similar to that of the T 3 topology given in Refs. [2, 13]. Equations (2.7) 

can be modified by defining e p = e p / sin (9 aand collecting the first and last terms of each 
equation. They become 



e~ T csc(e- T )(e T sin( e - T )P' T ) 



. 2r (sin^) 5 , 



sin6> 

e p '{Q 2 T + e- 2r Q 2 d } = 0, (3.1a) 

2T (sm9Q jd ) if ) 



e T csc(e T )(e T sin(e t )Q, t ),t — e 



sin 9 

2t tV 



+2(P; r Q, T - e- 2T P; d Q, ) = 0. (3.16) 
If we define a new time variable w by means of dw = e~ T csc(e _r )dT, or 



w = — In 



tan 



(3.2) 



then, using sin(e r ) = sech w, and multiplying through by e 2r sin 2 (e r ), Eqs. (3.1) reduce 
to 



( sin 9P'a) a OD / , » „ 

PL,,, - sech 2 w- - e 2P {Q% + sech 2 wQ 2 e } = 0, 



Q ww — sech w 



sin 9 
2 (sin 0) 



sin 9 



+ 2(P; W Q, W - sech 2 w P[ Q,e) = 0. 



(3.3a) 



(3.36) 



These equations are similar to those of the T 3 topology [2], with P^g and Q^q replaced by 
(sin 9P'q\qI sin 9, (sin9Q i e),e/ sin 9 and e _2r replaced by sech 2 w. If we look at (3.2), we 

see that near the big bang, w — > —00, and near the big crunch, w — > +00, sech 2 w — > e _2 l™L 
Given the similarity between Eqs. (3.3) and those of the T 3 case, it is not difficult to 
construct the Hamiltonian form of the action that gives them [14], 



d9 dw 



u p/ p; w + u Q Q ;W 
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+ sin # sech 2 «;[(P' e ) 2 + e 2P ' Q 2 g 



(3.4) 



Since sech 2 w — ■> e -2 !™! for u; — > ±00, one might expect that the solutions (2.14) would 
become AVTD at both the big bang and the big crunch. In this case the potential terms 
would become negligible and the Hamiltonian density would approach 



TT 2 

H= p ' 



2sin# 



+ 



e- 2P 'n 2 Q 

2sin# 



(3.5) 



which is, again, similar to the T 3 expression given in Refs. [2, 13, 14]. This reduced 
Hamiltonian has an exact solution of the form 



P' = Pq + In [cosh vw + cos ip sinh v w] , 
e p ° sinV'tanhww 



n. 



Q = Qo + 1 

1 + cos ^ sinh 

v sinh vw + v cos tp cosh vw 
cosh vw + cos -0 sinh vw 

Il„ = e p °vsmipsm6, 



sin6>, 



(3.6a) 
(3.66) 

(3.6c) 
(3.6d) 



where Po, v, Qo and ip are arbitrary functions of 9. 

If the solution (2.14) is to be AVTD at either the big bang or the big crunch, we have 
to rewrite P and Q in terms of w and see whether P' and Q take the form of P' and Q 
in (3.6). Since P and Q in (2.14) are functions of sin(e _r ) and cos(e _r ), and we have 
sin(e _T ) = sechw, cos(e _T ) = tanhty and P' = P — ln(sin#), we find that 



P' = In 



(1 - (3 2 ) sechw + (3 2 sin 2 9 cosh w 
a^l-13 2 sin 0(^1 -/3 2 tanhu; + l] 2 + f3 2 cos 2 9) 



(3.7a) 



-2a/3sin 2 0[^l -/3 2 tanhu; + 1] 
(l-/3 2 )sech 2 w + /3 2 sin 2 9 



(3.76) 



There are at least two ways to take the limit asw^ ±00; one is to assume that sech w — > 
and that 

f3 2 sin 9 cosh w 



P' -> In 



In 



a v / W? 2 (2 - /3 2 sin 2 + 2^/1 - /3 2 ) 
/3 2 sin 9 cosh tu 



a v /T^ 2 (2 - (3 2 sin 2 + 2y/l - (3 2 ) 



+ In coshty, 



Q 



2a, 



- — [1 + y/T=]Ptanhw], 



(3.8a) 



(3.86) 



which imply i; = 1, V = tt/2, P = ln^sinfl/av/T^p - /3 2 sin 2 + 2y/l - /3 2 }), Qo = 
—2a/ f3. However, the coefficient of tanhiu in (3.8b) is not e~ p °, which it should be from 



(3.6b). It is worse if we calculate the momenta. For example, Uq = e p Q w sin9, and 
using the expression (3.7b) we find that Uq calculated in this manner is not the same as 
(3.6d). 

While this is the simplest limit, it is possible to keep terms of order e~ w and still have 
a limit that has approximately the form of (3.6). Since sechw — > e - '™', coshw — > e + ' w ', 
and e ± l u 'l = cosh \w\ ± sinh \w\, and tanh m; — > ±1, the limit of P' as w — > ±oo is 



P' 



In 



ay/1 -/3 2 (2 - (3 2 sin 2 9 ± 2^/1 - (3 2 ) 
1 - f3 2 cos 2 e 



+ ln 



[l-/3 2 (l + sin 2 fl)] . 

cosh w =f 5-t sinh w 

1 — fj z cos 2 y 



(3.9) 



which implies cosip = =f[1 — /3 2 (1 + sin 2 6>)]/[l — f3 2 cos 2 0], u = 1, 
P = (1 - /3 2 cos 2 6»)/{a v /I^ 2 (2 - /3 2 sin 2 9 ± 2^1 -/? 2 )}. 
As w — > ±oo, tanhw — > ±(1 — ie -2 !™!) and Q becomes 

-2a/3sin 2 6'[± v /l -/3 2 (1 - ±e" 2 H) + 1] 



(l-/? 2 )e- 2 l»l +[3 2 sm 2 e 
and using, as above, e ±][w][ = cosh |«;| ± sinh \w\, we find that 



(3.10) 



Q 



-2a/?sin 2 0[(l + W1-/3 2 ) coshu; ± (1 + f^l - /3 2 sinh 



(1 - /3 2 cos 2 0) coshw ± [-1 + (3 2 (1 + sin 2 0)] sinhw 
-2af3sm 2 9(l + § ^1 - /3 2 ) 



1+1/2 v^Eg ± tanh 



(l-/3 2 cos 2 6>)[ 



^(l+sin^g)-! , 

1 =F 1-/32 cos 2 e tannw 



(3.11) 



which is of the form of (3.6b) with Q = -2a/3sin 2 0(1 + - (3 2 )/(l - (3 2 cos 2 0) and 
ifj as above. However, if we calculate Uq as before, it is not the same as (3.6d). 

The fact that these solutions do not become AVTD is, in fact, not surprising if we 
study the velocity v [2, 14], 



v = J-/u 2 pi +e- 2P 'U 2 
sin V ^ v 

that appears in (3.6). This quantity can be calculated for all values of w, and is 

. sechw 

v{9,w) = 7. 7, ; 

[(1 - (3 2 ) sech 2 w + (3 2 sin 2 9]{[ X /T^ft 2 tanh w + l] 2 + (3 2 cos 2 0) 

x ^{sinh«;[(l - (3 2 ) sech 2 w - (3 2 sin 2 9]([y/l - (3 2 tanhw + l] 2 + (3 2 cos 2 0) 
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+2sech«;[(l - (3 2 ) tanhw + ^1 - (3 2 ][(1 - (3 2 ) sech 2 w + (3 2 sin 2 6]} 2 

\ 1/2 

+4f3 2 sin 2 9{[y/l - /3 2 tanhw + l] 2 - (3 2 cos 2 6) 2 J . (3.12) 

As w — * ±oo, i> — » 1 and P' — > \w\. 

If we consider Eq. (3.3a), the solution will be AVTD if the terms involving 9- 
derivatives decay much more rapidly than those involving tu-derivatives. The term in- 
volving Q 2 9 , 

sech 2 «;e 2P 'Q 2 ,^e- 2 ^le 2 ^lQ 2 , 

does not necessarily decay rapidly asw^ ±oo. However, calculating we see that it 
goes to zero as w — > ±oo, so this term does not contribute in the limit. However, calculating 

e 2P Q 2 W and P' ww , we see that they decay as sech 2 w (which is consistent with the AVTD 
form from [3.6], where they should decay as sech 2 vw). However, (sin 9P' e )j/ sin6> becomes 
constant in w as w — > ±oo, so sech 2 w (sin 9P' e ) j g/ sin 9 is of the same order as P' ww 
and e 2P 'Q 2 w , and it is not surprising that the solution (2.14) does not become AVTD as 
w — > ±oo. Of course, if v were less than one, P' ww and e 2P ' Q 2 W would decay more slowly 
than sech 2 w (sm9P' e ) : g/ sin 9. 

As we mentioned above, the "spiky" behavior occurs far from the big bang and the big 
crunch. From Figs. 1 and 2 we see that the maximum deviation from a constant function 
of 9 occurs at r = — ln(7r/2), which is at w = 0. It is not difficult to show that the solution 
near w = is not AVTD, but a formal investigation of the behavior of v near 9 = 7r/2 and 
w = sheds some light on the origin of the spike at 9 = tt/2. 

If we consider v{9, 0), we find that 

2Vl-/?Wfl 
^° )= 1 + /3W0 (3 - 13) 



which varies from 2a/1 — (3 2 / (1 + (3 2 ) at 9 = 0, it to two at 9 = tt/2. The value at 9 = 0, it 
is always less than two, and for (3 near one the difference becomes large. Notice that the 
system becomes formally AVTD at w = 0, 9 = 0, in that both (P'g) 2 and (Q,o) 2 are 
zero at this point, while v = 2. If one could assume that this fact was relevant to the 
time development of the system, it would show that there is an indication that the spiky 
behavior of Refs. [2, 5], which is tied to v > 1 at specific points of the manifold, is present 
here also. However, in our case, "velocity dominated" at one point in 9 is meaningless, since 
the ^-derivative terms in (3.3) can be shown to be considerable at 9 = tt/2. Nevertheless, 
the existence of these velocity peaks may give us some clue about the true meaning of the 
spikes. 

One last property of our solutions is their behavior as functions of /?. As f3 — > 0, the 
Kerr solution (2.12) becomes Schwarzschild, and our Gowdy model becomes Kantowski- 
Sachs (compactified in the "radial" coordinate), and one of the singularities (w — > — oo) 
becomes a curvature singularity. For the equivalent of extreme Kerr, (3 = 1, both A and 
Q seem reasonable, although time independent, but P is singular. However, the inner and 
outer horizons become the same, and there is no Gowdy solution. Notice that t from (2.13) 
becomes a for all values of r. 
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IV. Conclusions and Suggestions for Further Research 

We have shown that the Kerr metric between the inner and outer horizons can be 
interpreted as an exact solution to the Gowdy models with S 1 x S 2 topology. The Gowdy 
solutions generated by this method are relatively flat functions of 9 at the horizons (the 
"big bang" at the inner horizon and the "big crunch" at the outer horizon) which grow a 
spike near 9 = n/2 in both the Gowdy functions P{9, r) and Q(6>, r). 

The singularities in the solutions are just lightlike surfaces where the solution passes 
from the coordinate patch suitable for the Gowdy interpretation to patches where we have 
true Kerr (black hole) metrics. One important point is that we have compactified in the 
5-direction which becomes the t-direction in the Kerr patches, giving solutions with closed 
timelike lines in the outer regions, a property shared with the Taub-NUT solution. Notice 
that if one compactifies the Kantowski-Sachs model (2.9) in the r-direction, there are also 
closed timelike lines in the Schwarzschild portion. 

The solutions we have found are not AVTD near the singularities, which can be shown 
directly from the solutions and the Hamiltonian form of the action which generates their 
equations of motion. The fact that the velocity that would appear in the AVTD solution 
is one seems to give a convincing argument why this is so. 

There are a number of directions in which this analysis can be extended. An obvious 
idea is to study numerical solutions of Eqs. (3.3), which might shed light on the structure 
of solutions for these models more general than the particular solution (2.14), and could 
give reasons for the form that (2.14) takes. 

Another possibility [16] would be to create new solutions to the models where the big 
bang and big crunch "singularities" become true curvature singularities. The technique 
for doing this is described in Ref [17]. 

One final way to extend our idea would be to consider other models which have 
cosmological sectors and "black hole" sectors separated by horizons. There are a number 
of metrics where one or the other of these sectors is well known, but where the other has 
never been interpreted as either a "black hole" solution or a cosmological model. The 
cosmological sector of some of these metrics can be compactified as we did with the metric 
(1.7) (or they may naturally admit a compact topology), leading to "black hole" sectors 
with closed timelike lines. The global geometry of some of these manifolds (including 
the Kantowski-Sachs model compactified in the r coordinate of Eq. [2.9]) has not been 
investigated in detail. 
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FIGURE CAPTIONS 



Figure la-g. These figures show the evolution of P — h, h = ln[a^/l — j3 2 sin(e~ r )], as 
a function of 9 for the solution given in Eqs. (2.14) for a = 1 and (3 = 1/2. Fig. la 
corresponds to r = —1.1447, Fig. lb to r = —1.144, Fig. lc to r = —1.1, Fig. Id to 
r = -0.75, Fig. le to r = - ln(yr/2), Fig. If to r = +5, and Fig. lg to r = +10. 

Figure 2a-d. These figures show the evolution of Q as a function of 9 for the solution given 
in Eqs. (2.14) for a = 1 and [5 = 1/2. Fig. 2a corresponds to r = —1.1439, Fig. 2b to 
r = -1.1, Fig. 2c to t = - ln(yr/2), and Fig. 2d to r = +5. 
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